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A theory for the relation between the spin polarization and luminescence in silicon is presented. 
The theory provides intuitive relations for phonon-assisted optical transitions between the conduc- 
tion and valence band edges. It is shown that an opposite behavior of longitudinal and transverse 
optical phonon-assisted transitions is responsible to recent experimental results of spin injection into 
silicon. The effects of spin-orbit coupling and doping are studied in detail. 
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Optical orientation and luminescence polarization are 
widespread techniques in studying the spin of electrons in 
direct gap semiconductors [H, Q. These techniques, how- 
ever, are not straightforward in indirect gap semiconduc- 
tors where optical transitions are mediated by electron- 
phonon interactions ■ In spite of a recent important 
progress in spin injection into silicon there is no 

theory that enables a quantitative analysis of the lumi- 
nescence from spin injected silicon. Such knowledge can 
be used to accurately determine the electrons' spin po- 
larization from the measured circular polarization (s- llj. 
Similarly, it can be used to infer the spin iiijection effi- 
ciency across ferromagnet / silicon interfaces |7H15| . 

In this letter we study the relation between the spin 
polarization of electrons in silicon and the circular po- 
larization degree of the luminescence. By invoking sym- 
metry arguments, we first provide concise selection rules 
for each of the phonon-assisted optical transitions in an 
unstrained bulk silicon. Then in the main part, these 
rules help in analyzing the polarized luminescence spec- 
trum calculated by a comprehensive rigid-ion model. Our 
theoretical results elucidate recent experiments of spin 
injection from iron to silicon 

The luminescence of free carriers in bulk silicon in- 
volves transitions between electrons from the six equiva- 
lent conduction band valleys along the A-symmetry axes 
and heavy or light holes at the top of the valence band 
(Fig. [Ha)). Due to the crystal translational symmetry 
and the indirect absorption edge, phonon emission or ab- 
sorption are required to offset a crystal momentum differ- 
ence of fco~0.85 X (27r/a) where a=5.43A is the lattice 
constant. The light intensity is proportional to. 
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and related matrix elements denote, respec- 

tively, radiation-matter and electron-phonon interac- 
tions where e is the light polarization vector and £ is 
the phonon mode: transverse-acoustic (TA), transverse- 
optical (TO), longitudinal-acoustic (LA) or longitudinal- 
optical (LO). The angular frequency of the photon is ujq 
and of the phonon is ug. The initial state at the con- 




FIG. 1: (a) Band structure and symmetry notations of sil- 
icon along the A axis. The dot-dashed and dashed arrows 
represent the two dominant contributions for phonon-assisted 
optical transitions. The lower left (upper right) inset magni- 
fies the splitting of the bands due to spin-orbit coupling close 
to the r25' (Tis) zone center region, (b) Phonon branches 
along the A axis and their corresponding symmetries, (c) 
Phonon dependent circular polarization properties for optical 
transitions of spin-up electrons from valleys whose axis lies 
along (cr^) or perpendicular {a'^) to the light propagation. 



duction band is given by \i) and the final state at the 
valence band by |/). The first and second terms in- 
volve, respectively, optical transition paths with inter- 
mediate states, whose crystal momentum is zero or 
kf). The dot-dashed and dashed arrows in Fig.[T](a) show 
the important respective paths where the vertical (non- 
vertical) arrows denote the photon (phonon) interaction. 
The circular polarization of the spin-dependent lumines- 
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FIG. 2: Transverse-phonon-assisted optical transitions of 
spin-up electrons from the x valley. The upper and lower 
diagrams correspond to different phonon modes whereas the 
left and right diagrams to circular polarization of x±iy. The 
numbers next to the arrows denote the relative amplitude of 
the corresponding transition. Note that if there are several 
paths to reach the same final state then the interference is 
constructive (numbers in brackets have the same sign). The 
relative contribution of each phonon and photon configura- 
tion is indicated at the top right corner of each diagram. The 
circular polarization ratio (crj : a'£) is 1:2. 



cence strongly depends on the phonon dispersion and 
symmetries along the A-axis. This information is shown 
in Fig. [IJb) . The TA and TO phonons share the same 
symmetry and therefore to first order their associated op- 
tical transitions are destined to have similar polarization. 
However, it will be shown that this symmetrical behav- 
ior breaks down due to the proximity of the LO and TO 
phonon branches around fco- 

Using the symmetries of the electronic states and of 
the phonon modes, one can derive concise polarization 
ratios between the left and right circularly polarized lu- 
minescence. Figure [IJc) lists the phonon-assisted cir- 
cular polarization ratios and the ensuing total circu- 
lar polarizations. These ratios correspond to transi- 
tions of spin-up electrons and they depend on whether 
electrons reside in valleys whose kg is parallel or per- 
pendicular to the light propagation axis (which is also 
the spin quantization axis). Figure [2] shows the cor- 
responding transverse-phonon-assisted optical transition 
diagrams from the x valley (perpendicular to light prop- 
agation) . Details about the construction of this and sim- 
ilar diagrams are found in the auxiliary material [l6| . 

The interference effects between the various intermedi- 
ate states is crucial for the analysis of the circular polar- 
ization. We focus on the resulting opposite polarization 
of the LO and the other phonon-assisted optical modes 
(Fig. IHc)). This suggests contributions of opposite sign 
to the circular polarization of emitted photons whose en- 
ergy is about Eg — huJop where Eg is the gap energy and 
ujop lulo ~ '-^TO (due to the nearly degenerate ener- 



gies of the LO and TO phonons). These photons will 
have smaller circular polarization compared with pho- 
tons whose energy is around Eg — Hlota in spite of the 
similar symmetry of TA and TO phonons. Other ef- 
fects of the interference are also instructive for analyz- 
ing the luminscence spectrum. First, the interference 
of zone center and fcg intermediate states is construc- 
tive (destructive) for transverse (longitudinal) phonons 
171 Il8| . This is the reason for the negligible intensity 
of the LA phonon-assisted optical transition. The sym- 
metry of LO phonons, however, precludes transitions to 
the Fis intermediate states [l9|. Consequently, the LO 
phonon-assisted optical transition and its circular polar- 
ization are dominated by intermediate states below the 
gap. The destructive interference with symmetry allowed 
intermediate states above the gap (e.g., the r2') is not 
sufhcient to destroy the LO phonon-assisted transition or 
to significantly alter its associated circular polarization. 

The robustness of the circular polarization and the rel- 
ative intensity of the various spectral lines are calculated 
using realistic modeling of Eq. ([Ij . The electronic states 
and energies are calculated via an empirical local pseu- 
dopotential method with spin-orbit coupling [20|]. The 
latter generates the correct energy splitting between T^ 
and r^" (Aso=44 meV). The phonon modes and their dis- 
persion are obtained by the adiabatic bond charge model 
mi. The calculated energy band structure and phonon 
dispersion curve along the A-axis are shown, respectively, 
in Figs, [ija) and (b). The radiation- matter matrix ele- 
ments are calculated using the electric dipole approxi- 
mation. The electron-phonon matrix elements are calcu- 
lated via the rigid- ion approximation 22] , 

gl.g2Xl,X2 

x{q-u^cos(Ag-x) + q-uisin(Ag-T)} (2) 

The wavefunctions are taken from the pseudopotential 
model and given by V3k,m = Eg, ,xj ^ei^Xj exp(z(gj+k)-r) 
where k, g, m and x denote, respectively, the wavevec- 
tor, reciprocal lattice wavevector, band index and spin 
state. Vq is the pseudopotential of wavevector q = 
Ak — Ag where Ak = k/ — k^ and Ag = gi — g2. 
yl=z-\//i/(2Mw£_Ak) where M is the mass of a silicon 
atom. The phonon displacement vectors, and ul, 
are the "in-phase" and "out-of-phase" motion of the two 
silicon atoms in the unit cell, respectively. These vectors 
are calculated via the adiabatic bond charge model. The 
equilibrium atom positions relative to the origin (mid- 
point) are given by ±t where r = (a, a, a)/8. The effect 
of the spin-orbit potential during a virtual transition is 
negligible and thus we consider only the local pseudopo- 
tential part in phonon-assisted transitions to intermedi- 
ate states. We have verified that the Elliot- Yafet spin-flip 
mechanisms and their interference are of minor impor- 
tance in our case of interest 23 , 24 1 . 



TABLE L Relative light intensities {Ir=Icr++Icr_) and circu- 
lar polarization degrees [P=(I^^—I„_)/ 1^) due to transitions 
with spin-up electrons. The light propagates along the +z di- 
rection. The light intensities are normalized with respect to 
the TO intensity from the x or y valleys. The total intensi- 
ties are the sum of contributions from all six valleys. The LA 
phonon results are not shown (negligible intensity). 



valley 
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total 


mode 


TO LO TA 


TO LO TA 


TO LO TA 


It 


1 0.115 0.086 


1.41 0.23 0.092 


6.82 0.92 0.53 


P[%] 


-32.3 5.3 -36 


0.01 50.1 0.7 


-18.8 27.7 -23.5 



Table U lists numerical results of relative light intensi- 
ties and their circular polarization due to transitions of 
spin-up electrons at the bottom of the conduction band 
to heavy and light hole states at the top of the valence 
band. The left (middle) column corresponds to the con- 
tribution from each of the four (two) conduction band 
valleys whose direction is perpendicular (parallel) to the 
light propagation direction. The relative total intensities 
from all six valleys (right column) are consistent with the 
luminescence spectrum in silicon (25| . The polarization 
values are consistent with those derived by symmetry ar- 
guments (Fig. [ijc)). The non-zero polarization of the 
LO mode from perpendicular valleys (~5%) is due to the 
small effect of transition paths via the V21 conduction 
band. In addition, the small deviation of the TA and 
TO total polarizations from 20% (as predicted by sym- 
metry arguments) generally comes from differences in the 
transition intensities of zone center and fco intermediate 
paths. The pseudopotential interpolation technique (V^ 
in Eq. ([2])), and the origin for the different intensities of 
TO and TA related optical transitions are spin indepen- 
dent and given in the auxiliary information [l6j . 

The spectral width of the circularly polarized region 
is strongly affected by doping. The previous results were 
derived for transitions between the extremal points of the 
conduction and valence bands. However, in a p-type (n- 
type) doped silicon one should consider all of the states in 
the valence (conduction) band that can take part in the 
transition as well as the ensuing larger set of intermediate 
states. We have integrated the transition probability of 
the various states in the Brillouin zone weighted by Fermi 
distribution at T =10 K. Fig. [3l[a)-(c) shows the result- 
ing spectra of Oj^ and ct_ for light propagation along the 
+z direction in 4x 10^^ cm^ n-type silicon, (a)-(c) denote, 
respectively, the contribution from the z, x and all six val- 
leys. The dominant peak just above 1065 meV originates 
predominately from the TO phonon-assisted transitions. 
Their associated circular polarization signal comes from 
the four transverse valleys (±a;, ±?/). About 15% of the 
dominant peak originates from LO associated transitions 
which are responsible for the small polarization of the 
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FIG. 3: Calculated polarized spectra in doped silicon at 
10 K. (a)-(c) show results of 4x10^* cm^ n-type silicon, 
(a) / (b) show the contribution from valleys whose axis is 
along/perpendicular to the light propagation direction, (c) 
shows the contribution from all six valleys, (d) shows the 
polarized spectra for 4x10^* cm'^ and 2x10^' cm"* p-type sil- 
icon. The former is magnified for clarity. The spectra in 
(a)-(d) consider the energy gap reduction with doping [25l]. 
(e) shows the polarization amplitude at the TA peak and TO 
dominated peak versus doping level (at the metallic regime). 



longitudinal valleys (±z). The opposite sign contribu- 
tions to the total circular polarization (Fig.[31[c)) results 
in an evident smaller polarization degree of the dominant 
peak compared with the less intense TA associated peak 
(around 1105 meV). Figure [3]Jd) shows the total spectra 
in cases of 4x10^^ cm'^ and 2x10^^ cm^^ p-type silicon. 
We notice two features when comparing the n-type and 
p-type spectra at a doping level of 4x 10"'^* cm'^ (Figs.[2Ic) 
& (d)). First, the p-type peaks are broader due to the 
position of the Fermi level; Ep^\^ meV below the va- 
lence band edge vs Ep^S meV above the multi-valley 
conduction band edge. Second, the p-type spectrum is 
less polarized due to stronger mixing of hole states away 
from the F-point. Both the broadening and reduced po- 
larization are evident in the spectrum of the 2x 10^^ cm'^ 
p-type silicon. Finally, we can identify a common rule 
shared by the spectra of all doping cases in the metallic 
regime (>2xl0^* cm'^). The polarization has its largest 
amplitude at the lower energy edge of a given phonon- 
assisted transition. At this photon energy, the initial and 
final states are from the extrema states whereas at higher 
energies the mixing of states plays a key role in lowering 
the polarization. 

The effect of spin-orbit coupling is most evident at 
heavily p-type silicon where the Fermi level is close to 
(or below) the split-off band. The proximity of the split- 
off band renders the holes state mixing effective already 
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~10 meV away from the F-point. Fig. [3je) shows the 
polarization values at the TA peak and TO dominated 
peak versus doping levels in the metallic regime. In n- 
type silicon, the polarization reduces only slightly with 
doping concentration. The initial polarization reduction 
of the TO dominated peak is due to the merging with 
the LO peak. This is essentially a broadening effect and 
it is effective when Ep>fiujTO-fvjJLO~'^ meV above the 
conduction band edge (>4xl0^^ cm^). In p-type cases, 
the peak polarizations decrease monotonically with dop- 
ing due to state mixing in the valence band. The latter is 
a general feature in semiconductors stating that emitting 
circularly polarized photons is not possible when Ep is 
well below the split-off band ^] . 

The results of Fig. [3] elucidate the measured polarized 
electroluminescence (EL) from spin injected silicon 0- 
[loj . The theory and experiment are consistent about 
the higher polarization of the TA related peak as well 
as for the polarization ratios between the TA and the 
TO dominated peaks. We find that the spin polariza- 
tion of injected electrons in Refs. is about 27% by 
comparing the measured and maximal achievable circular 
polarizations of the TA peak at 1105 meV. This relies on 
the ratio between the measured circular polarization of 
3.5% and the maximum attainable circular polarization 
of 13% at lO^^ cm^ p-type sihcon (Figjle)) [26i]. An- 
other point of consideration in Refs. ISMlOj is the <10 ns 



Auger lifetime of electrons in the p-type substrate [27 1. 
The EL measurements probe those electrons that recom- 
bine radiatively within the effective lifetime. Quantita- 
tively, about taItt of the electrons that reach the sub- 
strate experience radiative recombination where ta and 
Tr are, respectively, the effective (Auger dominated) and 
radiative lifetimes. The relatively short ta explains the 
mitigated spin relaxation effect in these measurements. 

In conclusion, we have studied the luminescence in 
spin polarized silicon. Useful circular polarization de- 
grees were provided for each of the phonon-assisted tran- 
sitions. The antipodal behavior of optical phonons was 
shown to be responsible for the differences in the circu- 
lar polarization of the luminescence peaks. The doping 
was shown to affect both the shape and polarization of 
the spectrum. Knowing the theoretical circular polariza- 
tion values of complete spin polarization and comparing 
them with measured values are imperative in determin- 
ing the spin polarization in silicon. These values are also 
instrumental in extracting the spin relaxation time or the 
spin injection efficiency across ferromagnet / silicon inter- 
faces. Finally, the theory sets a basis for future studies 
of the circular polarization in indirect gap semiconduc- 
tors due to strain (splitting the valleys and hole states), 
no-phonon and multi-phonon-assisted optical transitions. 
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